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Abstract 

We show that the delta function potential can be exploited along with perturbation theory to 
yield the result of certain infinite series. The idea is that any exactly soluble potential if coupled 
with a delta function potential remains exactly soluble. We use the strength of the delta function as 
an expansion parameter and express the second-order energy shift as an infinite sum in perturbation 
theory. The analytical solution is used to determine the second-order energy shift and hence the 
sum of an infinite series. By an appropriate choice of the unperturbed system, we can show the 
importance of the continuum in the energy shift of bound states. 
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I. INTRODUCTION 



Perturbation theoryi'^ii^i^'^i^'^i^i'i^ is a common tool in the study of eigenvalue problems. 
A classic eigenvalue problem is the time-independent Schrodinger equation, which can be 
written as 

Hi) = V^^ + V^ = E^, (1) 

2m 

where H is the Hamiltonian, m is the mass of the particle, the wavefunction is ip, and the 
eigenenergy is E. The problem is to find the allowed energies En and the corresponding 
wavefunctions ipn- The energies En are generally discrete and determined by the boundary 
conditions. Among the very few cases for which E and ip can be determined exactly are a 
particle in a box, a particle confined by a harmonic potential, and a particle confined by a 
Coulomb potential. The shape-invariant potentials for which exact solutions can be found 
are variations on these cases.— 

We denote the Hamiltonian for which an exact solution can be found by Hq and write an 
arbitrary Hamiltonian H as 

H = Ho + XH', (2) 

where A is a parameter and H' is the perturbing Hamiltonian. We denote the eigenvalues of 
Ho by -En*-"-* and the eigenf unctions by ipn \ and express the eigenvalues En of if in a power 
series expansion in A as = En^ + XEn^ + XEn^ + . . . with 

4" = mn'm, (3) 

p(o) p(o) • 

Equations ([3]) and (jlj) are valid when ipn^ is nondegenerate. From Eq. (jl]) we see that the 
expression for Ek^' is an infinite sum. If we choose the perturbation such that both H and 

(2) 

Hq are exactly soluble, we can obtain an exact expression for En which, when equated with 
the expression in Eq. (jlj), leads to an infinite series whose sum is known.— 

The (5-function potential can be used to fulfill this purpose. The solutions of Hip = Eip 
are the same as the solutions of Hqi/j = Eoip except where the ^-function is nonzero, and E 
can be found by matching the solutions across the 5-function. This procedure can provide 
some other insights as well. In Sec. II we discuss an example where the continuum states 
can change the sign of the contribution coming from the bound states. Another example is 



the effect of physical confinement on bound structures such as molecules. In this case we can 
study the effect of enclosing the 5-function(s) between two infinite walls which is explored 
in Problem 1 in Sec. VI. A similar extension of Problem 3, given in Sec. VI is also possible. 

In Sec. II we consider a box with infinite walls perturbed by an attractive 5-function 
potential placed anywhere within the box. Next we consider a box of finite depth so that 
the number of bound states is finite. In this case the impact of the continuum states on 

(2) 

the calculation of En becomes of primary interest. In Sec. Ill we use the simple harmonic 
oscillator with a (5-function potential at the origin to obtain an interesting infinite series and 
show how we can do the same for the one-dimensional hydrogen atom. The convergence 
of some other series is briefly considered in Sec. IV. We give a few suggested problems in 
Sec. VI. 



II. PARTICLE IN A BOX WITH DELTA FUNCTION(S) 

The 5-function potential has been used as a model for atoms, molecules, solids, liquids, 
the many body problem, and scattering.— The usefulness of the 5-function potential has 
been discussed by Lapidus^ in the context of stationary state perturbation theory. We 
place the potential at an arbitrary position, x = pL with < p < 1, in a box of length L, 
with < X < L. The general solution of Eq. ([1]) can be written as 

{A sin kx + B cos kx, (0 < x < pL) 
(5) 
C sin kx + D cos kx, {pL < x < L), 



where k = ^2mE /h?' and A, B, C, D are determined from the normalization and boundary 
conditions. Because ip{x) = at x = 0, we have B = 0. Hence we can write 

ip^x) = Asinkx, {0<x<pL). (6) 

We also have ip{x = L) = so that D = —C tan kL and 

^{x) = CsiYik{L-x), ipL<x<L) (7) 

where C = C / cos{kL). The continuity condition of the two functions at x = pL gives 

A sin ^(1 — p)L 



C sin kpL 
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(8) 



If we integrate Eq. ([T]) from — e to +e and take the limit as e ^ 

dx dx e^o J_ 



(9) 



It is clear that if V{x) is finite, ^ is continuous at x = 0. But when V{x) is infinite this 



agreement fails. In particular if V{x) = —X6{x) 

1 dip ^ 

doc d/Oc 



2mA , , , 



(10) 



Using this discontinuity condition at x = pL and Eq. ([8]), we obtain the eigenvalue 
condition, 



k sin kL 



2mA 
If 



sin kpL sin — p)L. 



If the perturbation is small, k can be expanded as 

k = + AA;(i) + A^A;^^) + . . . 



We write 



which gives 



= —e = — [fcw + \k^^) + x^k^^) + ...]', 

2m 2m 



From Eqs. ( ITTl) and ( |T2|) we obtain 



77(0) _ ^l{0)^ 

" "2m^ ' 

= ^2fc(°)fc(^, 
2m 

^2 



2m - 



fc«' + 2A;(o)A;(2) 



-11 



2m 



■sin2(A;(°)^^), 



A;(2) = + f!!^ [sin(A;(°VL) cos(A;(°V^)(l - 2p)] 



so 



A;{0) ' J,mf2 

that = n'^h^TT^/2mL^ and 

2 

= --sin^(n7rp). 



2m 



sin^(?27rp) [1 + 2nn(l — 2p) cot(?27rp)] . 



(11) 

(12) 
(13) 

(14a) 
(14b) 
(14c) 



(15a) 
(15b) 

(15c) 



(16a) 
(16b) 
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If we substitute the energy eigenvalues and eigenf unctions of the unperturbed Hamiltonian 
and write the perturbation as H' = —X6{x — pL) in Eq. (jl]), we find 

sin^(/7r]9) 



5m 2 
■ sm 



(nvrp) 



4 1^^?^ = sin^(n7rp) 



1 ^7r(l-2p) 
— + 2^ — cotinnp) 



n 



n 



From Eqs. ( ]16b[) and ( |T71) . we find 

rp 

If / and n are odd integers and p = 1/2, Eq. (ITSl) simplifies to 



P — -n? n? 



{I, n odd). 



(17) 



(18) 



(19) 



Equation ( fT9il is identical to Eq. (A3) of Ref. uA- To verify that the sum is correct and 
the series is convergent, we write Eq. f|T9l) as 



N 



P-V? 



— lim [ 

n N^oo I 



(2m + 1) - 72 ^ (2m + 1) + n. 



^ + lim 



_1 iifi Af-2±i Af+ii±i 

E ^- E ^ 



1-n 
2 



n + l 
2 



n + l 
2 



n Af^oo L ^ — ^ » — ^ qj ri'^ 



n + l P n + l ^" 

2 2 



(20a) 



(20b) 



(20c) 



because for finite n, — (n — l)/2 and + (n + l)/2 both approach as ^ 00. 

We can employ the first-order correction term to the wavefunction, which also has a 
sum-over-states form, to obtain a sum rule. We have 

A sin /cx, (0 < X < pL) 

,^^^li^L^^^KL - X) {pL<x<L) 



(21) 



with 



A' 



L sin kL cos 2kL{l — 2p) 
2 2k 



1 



(22) 



from the normalization condition. If we write k as in Eq. ([HD and write ^ = ^(°)+A^(l)+. . ., 
the first-order shift in the wavefunction can be found from 

sm X + Xk^^^'x cos X 



iIj^x) 



— sin k^^^'x 



cos k(o)L cos 2A;(o)L(l - 2p) 



1/2 



(23) 



1 + 



2A;(o) 



[ — 1)" cos4np7r 



+ X\ljk^^^xcosk^°^x + 0{\^) (24) 
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Thus the first-order shift ilJn\x) is 

If we use the fact that H'l^ = — sin rapvr sin /pvr, where = {tpl'^^lHlipn^), we can write 
Eq. §^ as 

n ^ , , A2' 2 sin npir sin /pvr sin ^ , ^ 



which leads to 



/2" 2 sin npvr sin Zpvr sin ^ [2^ nnx k^^^ , 
-\ / ^ — = A/ — sm -7-rv(-l) cos4n»7r 

+ y|A;«a;cos^, (27) 

with / and n odd and < x < pL. 

The effect of continuum states in second-order perturbation theory has been considered 
m Refs. llsIielHIlslQ. A one-dimensional attractive 5-function potential was considered in 



Ref. 



171 . This potential supports only one bound state^° and for a given perturbation, the 



contribution to En comes only from the continuum states. A one- dimensional square well 



18 



potential perturbed by another square well potential was discussed in Ref. 

In this context we discuss an example where the result of ii^n is a sum of a finite series 
for a finite number of bound states and an integral from the continuum. The continuum 
contribution is so large that it reverses the sign of the sum. We consider a finite square well 
potential, V = Q for —L < x < L and V = Vq outside for which there is a finite number of 
bound states determined by the magnitude of Vq.-'^^'^ In the following we will consider the 
bound state with the highest energy. The perturbation is an attractive 5-function at a; = 0. 
Because all the other bound states are lower in energy than the one we are considering, their 
contribution to the sum in Eq. (jlj) is positive. 

Equation ([T]) outside the well can be written as 

— ^(x) + ^ij{x) = 0, (28) 



with 



,-=^-^iVo-E). (29) 
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The corresponding solution is 



B exr)(—Kx) (x > 0), 
i^{x) = { (30) 
B exp{Kx) (x < 0). 



Equation ([T]) inside the well can be written as 



where 



The corresponding solutions are 



^ Mx) + k^tPix) = 0, (31) 
dx^ 



, r, 2mE , 
k' = ^. (32) 



A cos kx, even 
i^ix) = { (33) 

Asinkx, odd. 

The acceptability conditions for the wavefunctions are equivalent to requiring that the log- 
arithmic derivative of the wavefunction be continuous at x = ±L. For even states this 
condition leads to 



tanA;L = - = y^-l (34a) 



a=\l^^. (34b) 



If we plot a^L"^ — k'^L'^ and V fc^L^ tan^ V fc^L^ versus fc^L^, the intersections will give the 
scaled values of the bound state energies. Here all bound states are non-degenerate. Even 
and odd solutions alternate as the energy increases. The number of bound states is finite 
and equal to + 1 if Nn < 2aL < {N + l)n. Hence, by selecting a value of a, that is, Vq, 
we can control the number of bound states. 

If we perturb the system by placing an attractive 5-function potential at a; = 0, the 
wavefunction is given as 



ipi^x) 



A cos{kx) + B sin kx {0 < x < L) 

A cos{kx) - Bsmkx (0 > x > -L) (35) 
Cexp(— k|x|) (x > |L|) 
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The continuity condition at x = gives B = mXA/kh'^, and the wavefunction inside the 
well is 

^(x) = A\ coskx -\ gSi^^^l- (36) 

kh 

The continuity condition at x = +L gives 

k'^ t&nkL — kn = Xj-^lntenakL + k]. (37) 

If we substitute the values of k and k given in Eqs. (32) and (29) in Eq. (37), we obtain 
the desired result 

E tan ( V2mE^) - jE(Vo - E) = \jVo - E tan(V2mE^) + Ve] . (38) 

V nJ V2/i L h J 



tan ( V 2mEi°) ^) = , /^I, (39a) 



We expand E in powers of A, equate terms of same order of A on both sides and obtain 

Ei^) = ^ (39b) 



1 + 

Vo 1 



im3/2Le + ^^^^ + f^(2Ei°) - \/o) + ^(2^°^ - K,) 
^ V2h^ [h + LV2^ [1 + + fV2^ ' ^^^^^ 

where e = Vq — En^ . For the rath state to be bound En^ must be less than Vq. From the 
solution for En^ from Eq. (39a) we see that E^n^ <lVo, even for the most excited bound 
state. Given these constraints, we conclude that En < 0, which makes it clear that the 
integral from the continuum states has a greater contribution than the discrete sum. If we 
fix Vq and L such that there are only three bound states with energies Ei, E2, and E^, then 
wavefunctions corresponding to Ei and E^ are even and that corresponding to E2 is odd. 
The contribution from the bound states to the energy shift of E^ to second order is 

= (40) 

where Ai and A3 are defined in Eq. ( l35l) for the ground state and the second excited state 
respectively with a similar meaning for ki and k^. The relation between A and k is given by 

^^'^ = / ^ (41) 

/L±^sin(2A:i,3L) + cos2(%|^^) 



which is positive, where + and — signs are for Ai and A3 respectively. The contribution 
from the continuum makes the final answer as given by Eq. fl39cl) and is negative, showing 
the important role of the continuum. 



III. PARTICLE CONFINED BY A HARMONIC POTENTIAL 

We now consider the exactly soluble system of a simple harmonic oscillator perturbed by 
an attractive delta potential at the origin. The unperturbed Hamiltonian is = /2m + 
and the perturbation is H' = —\6{x). Equation ([1]) for the unperturbed system can 
be solved as 

(^2 / 1 z2\ 

—^{z) + U + - - - ^P{z) = 0, (42) 



rfz2-^' ' V 2 4^ 
where z = a/ 2'muj / hx and En = {n + l/2)huj, subject to the boundary condition that the 

wavefunction vanishes at ±00. The even and odd solutions of the unperturbed system are 

yi = e-^^/V,(| + i;i;^)even, (43a) 

y2 = ze-^-l\F, + 1 1 odd, (43b) 

where the true solution ip is a linear combination of yi and ?/2 (see Eq. (47)) and iFi is the 
confluent hypergeometric function (the Kummer function). Its power series expansion can 
be written as 

^ , , , a z a(a + 1) z"^ , , 

,f.(a;M) = l + ^^ + ^^ + ... (44) 

The solution can be written in the usual form as 

U (a, z) = cos ^ + ^) ^1 - sin ^{\ + ^) ^2, (45a) 

^^""'^^ = r(^-a) ['^^^(i + i)^i~^"'^(i + i)^2]' 

where Yi and Y2 are given by 



n = -^^^^yu (46a) 



We write the general solution for the perturbed problem as 



Ay, + By2 {x > 0) 
^ = <; (47) 
Cyi + Dy2 (x<0). 



These two forms of ip must match at x — 0, giving A — C. The wavefunction must be either 
symmetric or antisymmetric. The antisymmetric functions remain unchanged even after 
introducing the delta potential. The symmetric functions must change and are expressed as 



Ayi + By2, 
Ay I - By2. 



(48) 



The discontinuity condition of the (5-function at the origin gives B — —mX^A/k^, where 
^ = ^/h/2mOi), and the acceptable wavefunction takes the form 



iP^A 



yi 



1? 



-y2 



(49) 



The boundary condition ^0 — > as x — > oo leads to 



yi 



^2 = 0. 



(50) 



If we compare this form with Eq. (45a), we find that U{a,x) ^ as x ^ oo and thus we 
obtain 



V2 ^(i-t) tan(^ + ^) = -^. 
We express a and n in terms of En and expand the latter in powers of A and find 



(2n + -)hu;, 



moj r(n + \] 



0FV ^ r(n+l)' 



m 



r(n + |) 

r(n + l) 



V^(n + 1) - V'(n + -) 



(51) 

(52a) 
(52b) 

(52c) 



where 



ip{n + 1 



T'(n + 1) 



r(n + l)' 

1 rvti) 

2^ r(n + i)- 



(53a) 
(53b) 



On the other hand, if we substitute the values of the energy eigenvalues and eigenfunctions 

(2,) 

of the unperturbed Hamiltonian in Eq. (4) , we can express En as 

1 {2l)\{2n)\ 1 



" 27rh''^ 



22{n+o (n)2(^!)2 ^n-l) 



(54) 
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From Eqs. fl52cp and (IMI) . we find 



1 (2/)!(2n)! 1 


1 


rr(n + i)i 


2 p 






.r(n + 1)_ 





+ 1) - + -) 



(55) 



which is a new series with a known sum at the right hand side and may be of mathematical 
interest . 

We end this section with a brief look at the one-dimensional hydrogen atom. In this case 
we place the 5-function potential at x = a, so that the Hamiltonian is 

.2 j2 „2 



2m dx'^ \x\ 



(56) 



The eigenvalue equation is with E = —\E\ 



d'^i) 2a , I 2m ^. , , , 



dp"^ p 



h'\E\ 



(57) 



where a = m/2Ti \Eq\ = e^m/Ti k. The wavefunction is continuous aX p = ka and 



dip 


dip 




dp 


+ dp 





2m 
??\E 



■iP 



(58) 



at p = ka. For p < ka, ip is well behaved at the origin and vanishes at infinity. We have 



iP 



Here 



and 



Ape-PM{1 - a,2,2p), p < ka 
ip = Bpe-PU{l- a,2,2p) p > ka. 



M{a,b,z) = l + — +^-^ + ... 

(0)2 



U{1 - a, 2, 2p) = lim f/(l - a, 2 + e, 2p) 



TT 



lim 

e^o sin vre 



M(l - a, 2, 2p) 1 M(-a, -e, 2p) 



r(-a) 



2pr(l -a)r(-e)_ 



and {P)n = /?(/? + l)(/9 + 2) ■ ■ ■ (/3 + n + 1). We find the eigenvalue equation 



(59) 

(60) 

(61a) 
(61b) 



-M(l - a, 2, 2ka)U'{l - a, 2, 2ka) + M'(l - a, 2, 2A;a)[/(l - a, 2, 2ka) 

Tfl\ 

= ^M{l-a,2,2ka)U{l-a,2,2ka). (62) 
k 
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Now we can proceed to find the second-order energy shift. Although it is not instructive 
to show more details, outlining the procedure for the hydrogen atom is useful because it is 
a prototype for problems in higher dimension. In such situations, if we want to perturb the 
exactly soluble model by a 5-function, it is necessary to put the one dimensional 5-function 
at some r > rather than at r = 0, where it has no influence. 



IV. DISCUSSION 



To obtain a feeling for the convergence properties of the series that we have found, we 
present some numerical results. We first consider Eq. (fT9|) and rewrite it in the form 

Table [T] shows how the sum on the left approaches l/An^ in two sample cases, e.g. at n = 1 
and n = 5, with increasing number of terms. We can generate other series of this sort by 
taking other values of p{j^ ^) in Eq. (ITSl) . It is also interesting to determine what happens 
to the sum in Eq. fl63l) if the restrictions on I and n are removed. We can easily find that 
such a series converges to —l/An^: 

To verify, one may proceed as follows: 
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as, at the limit N ^oo, N±n^N. Consequently, the terms within parentheses tend 
to zero. 

We now turn our attention to Eq. (1271) which simplifies to the form 

EsinZpvrsin^ 1 nnx nx nnx , , 

— — = — - sm npn cos Anpn sm sm npn cos , (65) 
[V ~ n^) An^ L 2nL L 



where I and n are odd integers, < p < 1 and < x < pL. Unlike Eq. f[T8|) . there is an 
additional variable x so that we can find many interesting infinite series with closed form 
answers. For example, we consider p = 1/2, x = L/4, and = 1 in Eq. fl65l) that yields 

11111 1 71-2 , , 

+ TT-^ - ^r-^ - 1^ + + TTT^ + • • • = — — • (66) 



32-1 52-1 72-1 92-1 112-1 132-1 8 
If we group terms that are adjacent and have the same sign, we obtain 



00 ^ 



^ ^ 4kT2 ) ~ 

with no restriction on k. This series is interesting because it shows a sawtooth convergence. 
Alternate partial sums provide upper and lower bounds, which become gradually closer. For 
such series, it is well known that averaging two adjacent terms leads to better convergence. 
Table HI] shows some results for tt with a small number of terms in Eq. (1671) . We denote the 
sums by S{j) and construct S{j) = t^^^'^+f ^'^^^ to see how the average sum approaches better 
the desired value. By increasing the number of terms in Eq. fl67|) . we can obtain very good 
estimates of vr as shown in Table IIIII Note that the series (166|) or (!67|) is not commonly used 
to estimate vr. The well known Gregory-Leibniz series converges very slowly2^. Many other 



series for vr are available with varying convergence properties (see, for example, Ref. l2J). 
Series (67) may be a useful addition to such a list. 

We note finally that the infinite sums in Eqs. (|T8l) and (1671) for specific values of the 
parameters such as p and n or the variable x, have a common property: They have the same 
asymptotic behavior as the series f|T9|) . Hence, we can be sure about their convergence, in 
view of Eq. ([20]). 



V. CONCLUSION 



Exactly soluble potentials remain exactly soluble when we add a ^-function or a set of 
(5-functions (see Problem 3 in Sec. VI for an example). In the three cases we studied, the 
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problem can be solved in two different ways. An exact solution can be obtained by joining 
the piecewise exact results for the wavefunction at the ^-discontinuities. Or we can use 
perturbation theory in the strength of the 5-function to solve the problem to any desired 
order. Because the second-order perturbation theory result is usually in the form of an 
infinite series, we can use the exact solutions to determine the sums of such series, and 
exhibit their convergence properties. We can consider other exactly soluble potentials as 
well.— New infinite sums with known closed form answers may be found. Recently, an 
alternative scheme has been put forward using Green's function techniques.— 

By choosing appropriate examples, we can sometimes get insight into physics of certain 
problems that are not usually considered. The role of the continuum in the second order 
shift for the finite well is a case in point. In Problem 1 we explore the effect of a single 
attractive 5-function. We have seen that below a critical length, the bound state of the 
^-function disappears. Two attractive (5-functions can be used to set up a toy model for the 
hydrogen molecular ion. Placing walls in this case on the two sides allows us to study the 
effect of confinement on the binding of a molecule. We can imagine carrying this further 
and placing a periodic array of 5-functions interrupted on each side by a wall. The resulting 
violation of Bloch's theorem would make an interesting study. 

VI. SUGGESTED PROBLEMS 

Problem 1. In Sec. II we considered only the positive eigenenergies for the (5-function 
potential in a box. The attractive 5-function by itself will have a negative energy bound 
state. In this problem we explore what happens when the 5-function is placed in a box. 
Consider a 5-function at the center of a box of width 2L. Repeat the calculation of Sec. II, 
but for a negative energy bound state and show that for a 5-function of strength A, the 
negative energy bound state will disappear if L is smaller than Lc = Ti' jraX. 

Problem 2. Consider the perturbation of a particle confined in a one-dimensional box of 
length L by an attractive 5-function potentials at x = pL with < p < 1. The wavefunction 
is given as 



where A is a normalization constant. Explore the wavefunction in the region pL < x < L 




(68) 
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to obtain the first order correction to wavefunction. Also derive the same from perturbation 
theory and show that this leads to an interesting sum rule. Is the series convergent? 

Problem 3. Consider a particle confined to a one-dimensional box of length L perturbed 
by two attractive 5-function potentials at x = pL and x = qL with < j9 < 1 and < g < 1. 
(a) Find the second order energy shift and show that it is always negative when p = \ 
and q = J- (b) Do the same using Rayleigh-Schrodinger perturbation theory and compare 
the result with the above. Does this series converge? (c) Find the first order correction to 
the wavefunction and show how this correction gives three important sum rules which are 
convergent. 

Problem 4- Consider a simple harmonic oscillator and perturb it by an attractive delta 
potential at the origin. Use perturbation theory to find the first order correction to the 
wavefunction and show that this correction leads to another sum rule. 
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Tables 



No. of terms 


Results for n = 1 


Results for n = 5 


10 


0.2273 


-0.013 


50 


0.2451 


0.0051 


100 


0.2475 


0.0075 


1000 


0.2498 


0.00975 


10 000 


0.249975 


0.009975 


100 000 


0.2499975 


0.0099975 



TABLE I: Convergence study of Eq. (63) in two sample cases 



No. of terms 


Results for tt 


Average estimate 


[j] 


[sm 


[S{j)] 


1 


3.33 




2 


3.07 


3.20 


3 


3.18 




4 


3.12 


3.149 


5 


3.16 




6 


3.13 


3.144 


7 


3.15 




8 


3.13 


3.143 


9 


3.147 




10 


3.137 


3.142 


19 


3.1429 




20 


3.1404 


3.1417 



TABLE IL Convergence study of Eq. (67) for estimates of tt 
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j 


S{j) 


10 


3.142 


100 


3.141 593 


1000 


3.141 592 654 


10000 


3.141 592 653 59 


100 000 


3.141 592 653 589 793 8 



TABLE III: Gradually improved estimates of tt from Eq. (67) 
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